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ABSTRACT 

The ice impact problem on the engine blade made of layered composite is simulated. 
The ice piece is modeled as an equivalent spherical object and has the velocity opposite to 
that of the aircraft with direction parallel to the engine axis. Near the impact region and 
along the leading edge, the blade is assumed to be fully stressed and undergoes large 
deflection. A specified portion of the blade around the impact region is modeled. The 
effect of ice size and velocity on the average leading edge strain are investigated for a 
modified SR-2 model unswept composite propfan blade. Parametric studies are 
performed to study the response due to ice impact at various locations along the span. 
Also, the effects of engine speed on the strain and impact displacements are discussed. It 
is found that for a given engine speed, a critical ice speed exists that corresponds to the 
maximum strain and this critical speed increases with increase in the engine speed. 
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LIST OF SYMBOLS 


A, B coefficients in the displacement recursive relation 
A',B' coefficients in the velocity recursive relation 
a , b lower and upper radii for the impact modeling region 

d chordwise nodal distance from the leading edge 
Ell/ E 22 Modulus of Elasticity in 11 and 22 directions (psi) 

F, G displacement responses to unit initial displacement and unit initial velocity 

F', G' velocity responses to unit initial displacement and unit initial velocity 

/ impact force 

fi impact force on i th nodal 

f av impact force normal to blade chord (lb) 

Gj2 Modulus of Rigidity in 12 direction (psi) 

h time step (sec) 

[K] element stiffness 

Ki i th component of generalized stiffness 

i th component of generalized mass 
n number of materials 

P, i th component of generalized force 

r equivalent radius of ice piece 

T impact load duration (sec) 

t, T time (sec) 

t e element thickness (inches) 

t n time at n th time step (sec) 

V rel impact relative velocity (inckes/sec) 

\ element nodal coordinates (Figure 4 and Eq. 1) 
x i >yi ( or 

1 nodal coordinates with respect to impact node (Figure 6a) 
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Vj2 Poisson's ratio in 22 direction 

0 impact angle relative to the blade chord 

p density ( Ib.sec^jin 4 ) 

(Tyieid element yield stress 

coi natural frequency (radians) 

%i i th modal coordinate 

£,• velocity of i th modal coordinate 

acceleration of i ^ modal coordinate 
£ i n i th modal coordinate value at n th time step 
%i >n velocity value of i * modal coordinate at n th time step 
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INTRODUCTION 


At high altitudes, when aircraft flies through clouds of super-cooled water droplets, ice 
formation occurs on forward facing structural components. One such component is the 
engine inlet. With time, the ice accretes on the inlet and eventually sheds due to 
structural vibrations. A schematic of this phenomenon is shown in Figure 1. As a result, 
blocks of ice travelling at the speed of aircraft impacts the engine blades rotating at high 
RPM. This process may cause severe damage to the blade and subsequently to the engine. 
In order for the blade to sustain the ice impact, it is necessary to properly account for 
these constraints during the design. 

Fibrous composites are ideal for structural applications such as high performance aircraft 
engine blades where high strength-to-weight and stiffness-to-weight ratios are required. 
These factors along with the flexibility to select the composite layup and to favorably 
orient fiber directions help limit the impact damage and stresses arising from large 
rotational speeds. 

The objective of this paper is to simulate ice impact on an engine blade. SR-2 unswept 
composite propfan blade is considered for local damage analysis. The impact analysis is 
carried out by modifying the foreign object damage option provided in Reference [lj. It is 
assumed that the damage is severe in the local region and hence only a specified portion 
of the blade around the impact region is modeled. Furthermore, large scale deflections 
accompany the impact event in the impact zone [2]. In order to simulate this behavior 
near the impact region and along the leading edge, the blade is assumed to be fully 
stressed and undergoes large deflection. This is accomplished in the code by modifying 
the membrane stiffness in the radial direction to reflect a fully yielded condition. A 
detailed justification for these assumptions is also given in Reference [2]. Parametric 


4 



studies are performed to study the response of the blade due to ice impact at various 
locations along the span. The effect of ice size and velocity on the average leading strain 
and maximum impact displacements are investigated. 

IMPACT ANALYSIS 


Geometry of Ice Impact 

The geometry of ice impact on the leading edge of the blade is shown in Figures 2(a,b). 
Impact velocity direction relative to the blade is a function of aircraft speed and the 
rotational speed of the blade (Figure 2a). The resulting impact force not only depends on 
the magnitude and direction of relative velocity but also on the mass of ice piece. 
Depending on the spacing (i.e., number) of blades and relative velocity as shown in 
Figure 2b, only a portion of the ice piece hits the blade . The ice piece that is approaching 
the blade under consideration is sheared off by the adjacent blade and only a part of it 
impacts the leading edge. Effectively, the size of the ice piece that finally impacts the 
blade depends on the blade spacing, blade speed and aircraft speed. 

Modeling of the Blade 

The damage due to impact is considered to be highly localized and hence only a local 
patch: portion of the blade around the impact region (i. e., a specified portion along the 
span and half of the blade along the chord as shown in Figure 3) is modeled. In the code, 
the blade geometry is input in the form of finite element grid and nodal thicknesses. The 
spanwise impact region is specified with two parameters, namely, lower and upper 
bounds of radial fractions, a and b. However, the modeled impact span is approximated 
to be the region between the two finite element radial stations that are defined in the 
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blade geometry and are nearest to a and b. Along the chord, only half of the blade is 
considered for impact analysis. The finite element used is similar to NASTRAN (TRIA3) 
three node triangular plate element [3,4]. A total of 35 nodes and 48 elements are used 
(Figure 3). The impact is considered at the midpoint of the local patch (finite element 
node 16) along the leading edge. All the edges except the leading edge are considered 
fixed as the impact response beyond these boundaries is negligible. 

Large scale deflections usually accompany the impact event. This phenomenon is 
simulated in the elements that are close to the impact node by modifying the membrane 
stiffness in the radial direction (to reflect fully yielded condition) and zeroing-out the 
spanwise bending stiffness. These special elements are shown as shaded in the Figure 3. If 
CT yield is the effective yield stress of the element in the radial direction, w\, wi and are 
out of plane displacements at 3 nodes, the modified stiffness associated with these 
degrees of freedom is given by [2], 


[K] = |^yieid t e 

where t e is the element thickness and r 2 and y 3 are the element local coordinates (Figure 
4). The terms of this stiffness matrix replace the spanwise bending stiffnesses of linear 
elastic triangle. The stiffness in Equation (1) reflects the perfectly plastic condition with 
constant stress (<7 yield) in the radial direction. 

Figure 5 shows the material configuration of the composite blade. In each finite element, 
the skin (material I) is always present. The remaining materials exist only if the blade 
thickness permits them. Materials (n-1), (n-2) are dropped from the layup of an 
element in that order. If the element thickness can accommodate all the materials then 


1 -1 0 
-1 1 0 
0 0 0 


( 1 ) 
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the core (« th material with variable thickness) is used to fill the remaining excess 
thickness. 

Modeling of Ice Impact Loading 

The impact force on the blade arises from the momentum of the ice piece. The temporal 
and spatial distribution of this impact force is carried out with the following 
considerations [2]: 

a) Only the impact force component normal to the blade chord will cause the local 
deformation. 

b) Load duration, T, of impact is the time taken by the equivalent spherical ice piece 
(radius, r and relative velocity, Vrel) to squash-up, i. e., T= 2 r/V Te j. 

c) Impact force is distributed using a parabolic distribution function over all finite 
element nodes within one diameter (equivalent) of the ice piece from the impact 
node (see Figures 6a and 6b). 

_ fav t* ~ y»/2r ] ( 2 ) 

1 Z [l-x?/4r2][l-y/2r] 

i e n 

where /,• is the nodal force, fav is the component of the impact force normal to the 
blade chord, x; and y; are the local coordinates with impact node as the origin (Figure 
6a). 
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d) The impact wave travels along the chord at the same velocity as that of ice impact. As 
a result, i th nodal loading start time, T,- , is defined by the ratio of nodal chordwise 
distance, d from the leading edge and the ice relative velocity, V re i, i.e., T; = d/V re i, 
(see Figure 6a). 

e) Blade flexibility reduces the instantaneous ice impact loads due to reduction in the 
relative velocity. However, deflection of the blade results in ingestion of extra ice and 
therefore, there is an increase in the total momentum exchange. These factors are 
accounted and nodal forces are adjusted accordingly (see Reference [2] for details). 

Transient Response Analysis 

Transient response of the local patch is obtained by the process of modal integration. The 
modal analysis is based on the principle that statically deformed structure can be 
described by a combination of that structure's mode shapes. Once the coefficient of each 
mode's participation is calculated, the impact event can be simulated by the integration 
of a series of linear steps through time. The deflected shape due to ice impact is basic; 
therefore, only the first five modes of the local patch are used to compute the response. 
Also, damping is assumed to be zero as very little structural damping occurs during the 
impact event. 

Upon diagonalizing the governing equations of motion of undamped system [5], the 
uncoupled equations with each modal coordinate, £»', satisfies the second order 
differential equation. 


& + u? Zi - - V^W) 

Mi 


(3) 
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where a ),- 2 = M,- and K,- are the i th diagonal elements of the generalized mass and 

stiffness matrices, and Pi<t) is the i th component of the generalized force at time t. 

The general solution of the above equation, expressed in terms of arbitrary initial 
conditions, £> and £> at t = f„, and a convolution integral of the applied load, is, 

5<0 = KM«) 6 .» + G(t-r)Pi(t)dx ( 4 ) 


where F and G, respectively, the solutions with unit displacement and unit velocity 
initial conditions. 

Assuming that load varies linearly between t n and t n+ j , solution at time t n+ j can be 
obtained in closed form as (see Reference [6]), 


5,n+i = F(ft)5> + + A Pi >n + B P,>+i (5a) 

5>+l = F'(ft)5> + G'(fc)£,„ + A' P,> + B 1 P I>+1 (5b) 


where h = f n+I - 1„ , F(/i) = coscof i, C(7i) = (1/©,-) sin^/i, F'(/i) * sincu/i and G’(h) = 
coscOjh for undamped system [5]. 


For zero initial conditions (<5o = = 0) , the coefficients A, A’, B and B are given by 


A = 


1 

sin a)ih 

AW 

[ o)ih 


- cos 0 )ih\,B 


1 j x- sin_M 


(6a) 
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(6b) 




cos a )jh 

(Ojh 




The response at all time steps can be obtained by applying equations for &,n + 1 and £>+ 1 
recursively. At each time step, the modal coordinates, can be transformed to nodal 
displacements using modal matrix. From these displacements, large radial strains, £^, of 
the leading edge elements are computed with the relation [2], 



The entire calculation proceeds very quickly, due to the uncoupled nature of the modal 
responses. Among the strains at all the time steps, the maximum is picked as a 
representative quantity for the local impact damage. 


RESULTS AND DISCUSSION 

The SR-2 model unswept propfan blade [7] is chosen to study the ice impact and the 
planform of the blade is shown in Figure 7. The setting angle (orientation of the blade 
chord with respect to the plane of rotation at 75% span) is 57° and the number of blades 
are 8. In the present study, the composite material layup: Titanium skin/Graphite-Epoxy 
(±45)/Titanium core is used (Figure 8). For these materials, the properties are given in 
Table 1. The effective yield stress for element stiffness modification is assumed to be 206 
ksi. 
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Effect of Ice Size and Ice Velocity 


The impact in the region 50-90% of the span is known to cause severe damage to the 
blade [2]. Hence, lower and upper limits of the modeling region are taken as 50% and 
90% blade span for the impact analysis. For this modeling region, the impact radius is at 
9.55" (73.5% span). These data are presented for three engine speeds, 3000, 5000 and 8000 
RPM in the form of contour plots (Figures 9a-c). For all the engine speeds, the maximum 
strain corresponds to r = 0.8". However, the corresponding critical ice speed varies 
depending on the engine RPM. As an example, at 3000, 5000 and 8000 RPM the peaks are 
respectively at 100, 130 and 190 knots. The associated peak values of the strains are 1.66%, 
7.61% and 14.22%. From these plots, a safe range of ice impact parameter, namely, ice 
speed and ice size that can cause a specified local damage to the blade can be determined. 
Once the designer determines the tolerable strain along the leading edge of the blade, the 
combination of ice size and speed are defined by the region to the left of the level curve 
labelled with that strain. 

For the maximum ice size, a typical variation of strain with ice speed is shown in Figure 
10. The strain is zero at ice speeds 0 and 211 knots. But it reaches a maximum value at 
100 knots. This is due to the fact that only impact force normal to the chord (/ sin 6 ) gives 
rise to this strain. When ice speed= 0, the impact force if) is zero and when ice speed = 
211 knots, the impact angle ( 6 ) is zero. Initially the value of 6 is equal to the stagger 
angle (55° at the impact radius) and gradually decreases to 0° as the ice speed increases 
(Figure 10). 

Figure 11 shows the the effect of impact radius on the strain. The strain is nearly zero 
when the impact takes place below 53% span (approx. 7.5" radius) but steeply increases as 
the ice impacts towards the blade tip. This is expected as the blade is relatively thin closer 
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to the tip. Also, the blade velocity, consequently, impact force increases as the impact 
radius increases. 

Effect of Engine Speed 

With engine RPM as the parameter (1000-10,000 RPM), strain plot shown in Figure 10 is 
repeated in Figure 12. It can be seen that the variation of strain is smooth for engine 
speeds up to 4000 and fluctuates for the higher speeds. Flowever, the average trend 
remains the same for all RPM and the peak strain is larger for higher engine speed. The 
probable cause of these local fluctuations is due to discrete nature of the nodal force 
distribution, i.e., the finite element nodes which are in the force distribution boundary 
are suddenly loaded when the impact wave reaches it (Figure 6a). 

Relative magnitudes of typical impact (out-of-plane) displacement distributions for three 
engine speeds (3000, 5000 and 8000 RPM) are shown in Figure 13. The data corresponds to 
Ice speed of 200 knots and ice radius of 0.8". The displacements along the leading edge 
remain same and gradually decreases to zero toward the center of the blade. The 
maximum displacement for 5000 RPM is higher than that for 8000 RPM as the analysis 
ice speed (200) does not correspond to the critical one in both the cases. 


CONCLUSIONS 

The ice impact on an aircraft engine blade made of layered composite is simulated. From 
the results presented in this paper, the following conclusions can be drawn: 

1. The largest ice piece results in maximum average leading edge strain. 
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2. For a given engine speed, a critical ice speed exists that corresponds to the maximum 
strain. This critical speed increases with increase in the engine speed. 

3. The leading edge strain increases steeply with increase in impact radius. 

4. The strains generated from this analysis can be effectively used for structural 
tailoring of engine blades. 
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Table 1: Properties of the Blade Constituent Materials 



i 

Figure 1. A Schematic of Ice Impact on an Engine Blade 
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Figure 2. Geometry of Ice Impact on an Engine Blade 
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Figure 4. Finite Element Local Coordinate System 



Figure 5. Symmetric Material Layer Configuration 
in the Composite Blade 
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Figure 6b. Parabolic Force Distribution in the Region n 
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Figure 7. A Schematic of the Planform of SR-2 Model Propfan Blade 



Figure 8. Material Layup of the Unswept Blade Example 
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Figure 9. Average Leading Edge Strain Contours as a Function of Ice Speed and Ice Size 
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Figure 10. Variation of Average Leading Edge Strain and Impact Angle with Ice Speed 
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Figure 11. Effect of Impact Location on Average Leading Edge Strain 
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Figure 12. Effect of Ice Speed on Average Leading Edge Strain with RPM as Parameter 
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Figure 13. Absolute Out-of-Plane Displacements of Impact Region, Ice Size = 0.8",Ice Speed = 200 knots. 
Impact Location: 73.5% Span 
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